Condensate-free superfluid induced by frustrated proximity effect 
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Since the discovery of superfluidity in ^He and Landau's phenomenological theory, the relationship between 
Bose condensation and superfluidity has been intensely debated. ^He is known by now to be both superfluid and 
condensed at low temperature, and more generally, in dimension D>2, all superfluid bosonic models realized in 
experiments are condensed in their ground- state, the most recent example being provided by ultracold bosonic 
atoms trapped in an optical lattice. In this paper, it is shown that a 2D gas of bosons which is not condensed at 
T = can be achieved by populating a layer through di frustrated proximity effect from a superfluid reservoir. 
This condensate-free bosonic fluid is further shown to be a superfluid with incommensurate correlations. 



The low energy properties of many condensed-matter sys- 
tems are most naturally described in terms of bosonic par- 
ticles [1, 2]. This is of course the case of ^He [3, 4] and 
of cold bosonic atoms loaded in optical lattices [5], but ef- 
fective bosonic models have also provided an accurate de- 
scription of several systems such as superconducting thin 
films [6], Josephson junction arrays [7], or quantum magnets 
in a field [8]. At low enough temperature, most bosonic sys- 
tems exhibit two remarkable quantum effects: Bose-Einstein 
condensation (BEC) and superfluidity. These two manifesta- 
tions of quantum coherence often come together, but this need 
not be the case: for instance, free bosons condense but are not 
superfluid, while in principle superfluidity does not require a 
condensate, as first pointed out by Landau in his phenomeno- 
logical theory of superfluid ^He [9]. However, to destroy the 
condensate while keeping a system superfluid at T = in re- 
alistic situations turns out to be tricky. It is well-known [10] 
that interactions induce quantum fluctuations and generically 
tend to deplete the condensate in the ground- state (GS), but 
they usually do not empty it completely. For instance, in "^He, 
despite the very strong short-range repulsion, a condensate 
fraction of po ^ 7% of the total density remains [11, 12]. 
Generally speaking, interactions are known to deplete com- 
pletely the T = condensate only for ID geometries be- 
cause of diverging fluctuations and for long-ranged interac- 
tions in the 2D Bose Coulomb liquid model with ln(r) inter- 
actions [13], a model with no experimental implementation 
so far. Increasing further interactions will ultimately induce a 
crystallized phase, which is of course neither condensed nor 
superfluid. A realistic example for D> 2 where superfluid- 
ity could occur independently of BEC in the GS of a bosonic 
system is clearly missing. Here we explore a new route, using 
geometrical frustration to achieve a realistic condensate-free 
superfluid in D= 2. In the context of quantum magnetism, 
frustration has long been identified as a possible mechanism 
to produce spin liquid phases [14, 15], but it has only been 
realized recently that for bosonic models it could also induce 
new states of matter such as the supersolid state [16, 17] or 
the putative Bose metal phase [18, 19]. 

In this Letter, we show that the GS of a 2D bosonic fluid 
created by a frustrated tunneling from a reservoir (as depicted 
in Fig. 1) can be superfluid and not condensed. This conclu- 




FIG. 1 : Schematic picture for the frustrated bilayer made of two non- 
equivalent planes A and B which are coupled by a frustrated hopping 
term t±. The thick arrows depict the spin- 1/2 degrees of freedom 
[equivalent to hard-core bosons, see the text], which have their trans- 
verse component AF ordered [bosons condensed at (tt, tt)] in A layer 
while geometric frustration prevent AF order in B layer which, at the 
classical level, remain fully polarised [insulator for the bosons]. 

sion is quite general and only relies on the frustrated character 
of the coupling. It thus applies in principle to a variety of 
systems such as optical lattices loaded with bosonic atoms or 
frustrated dimers in a strong magnetic field like BaCuSi206, 
as we will discuss at the end of the paper. The starting point 
is a simple 2D model of hard-core bosons on a square lattice: 

n = t ^(ajar+^ii + h.c.) - /i ^ nr, (1) 

r,r|| r 

where ry runs over the two basis vectors of the square lat- 
tice. The hard-core constraint is essential to get superfluid- 
ity because it yields a linear excitation spectrum, which im- 
plies a finite critical velocity [20]. On the other hand, such 
an infinite local repulsion is expected to deplete the conden- 
sate with respect to free bosons. To quantify these effects, 
it is expedient to map the problem onto a spin model using 
the Matsubara-Matsuda transformation [3]: = nr — 1/2, 
= al, S~ = ttr, which, up to a constant, leads to an 
antiferromagnetic (AF) XY model in a transverse field. 

r,r|| r 
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with J = 2t and H = ji. The basic physics already appears at 
the classical level. For large negative /i, the system is empty 
(all the spins point down: {S^) = —1/2). At /i = /ic = — 4t 
(H = He = —2 J), the bottom of the quadratic single particle 
dispersion band £:(k) = 2t (cos kx + cos ky) — ji vanishes at 
ko = (7r,7r). This momentum defines the condensate mode 
in which bosons start to accumulate. The condensate density, 
defined on a lattice with N sites by 

PO ^ ^(4o«ko) = ^ E«^'°^"""^'^(^r+5r-), (3) 
r,r' 

is, in the magnetic language, equal to the AF order parameter 
in the plane perpendicular to the field. The superfluid density 
ps is proportional to the stiffness (or helicity modulus) T 
of the system, i.e. the second derivative of the energy with 
respect to a twist angle enforced at the boundaries [21]. At the 
classical level, the AF coupling between the XY components 
of neighboring spins induces a planar long-range Neel order, 
which in the bosonic language leads to non-zero condensate 
and superfluid densities given by po = Ps = — p)^ where 
p = 1/2 + /i/8t is the bosonic density [20]. So, at this level 
of approximation, the condensate and superfluid densities 
are strictly equal. Quantum fluctuations slightly reduce the 
condensate and increase the stiffness [20], but in the simple 
model of Eq.(l), the condensate persists at all densities. So, 
a single layer of hard-core bosons is a good prototype of 
a bosonic system in D> 2, with both a condensate and a 
superfluid density at T = 0. 

Bosonic bilayer — We want to populate a layer through 
quantum tunneling from a superfluid reservoir. Therefore, we 
consider two copies A and B of the model (1) 

H-^ = t ^{ala,+r, + h.c.) - ^i^nf (4) 

r,r|| r 

= t Y.{hlK+r, + h-c.) - (Ai - A) ^ , (5) 

r,r|| r 

where A > is an energy barrier. The critical chemical po- 
tentials now talce different values: on the A layer, pc = — 4t, 
whereas on the B layer, pc = — 4t+ A, and if < /i—/ic < 
the A layer has a finite density while the B layer is empty. Let 
us now consider, in such a regime, the effect of a frustrated 
coupling between the layers 

H-^" = t^Y.(^l+rA + bla,^,J, (6) 

r,r_L 

where r± runs over the four vectors coupling one site of the 
B layer to its four nearest neighbours in the A layer. As de- 
picted in Fig. 1, in the magnetic representation it is easy to see 
that such a transverse frustrated coupling leads, at the classical 
level, to a vanishing local xy field at each site of the B layer, 
and therefore prevents bosons from A to tunnel into B. Such 
a classical decoupling has to be contrasted with the case of a 
direct unfrustrated tunneling of the form t± Y2ri^lK + ^la^,) 



for which a local field is immediately induced in the B layer, 
leading to a finite density p^. 

Coming back to the frustrated case of Eq. (6), since the 
local field vanishes, the A layer has no influence at the 
classical level on the B layer, which remains empty. This 
cannot be true however when many -body effects and quantum 
fluctuations are included. Indeed, in the bosonic language, 
a wave-function with particles only in the A layer is not 
an eigenstate when one includes the interlayer hopping of 
Eq. 6. So particles have to be present in the B layer. This 
does not mean however that there is a condensate. Indeed, 
according to field theory, a direct (linear) coupling between 
condensate order parameters on A and B is incompatible with 
the symmetry of the lattice [22]. In our case of non-equivalent 
layers, this implies that a condensate in the B layer will 
not develop immediately when a condensate appears in the 
A layer. As we shall see, treating quantum fluctuations at 
the level of linear spin-wave (LSW) theory indeed leads 
to a small bosonic population in the B layer which is not 
condensed. We now explore the very peculiar properties of 
this bosonic gas. 

Linear spin-wave corrections — To treat quantum fluctu- 
ations, we start from the spin representation and perform a 
Holstein-Primalcoff transformation [23] after a local rotation 
of the spins. In the classical GS, the spins in the A layer make 
an angle ^Oa = ± arccos(— /i//ic) with the z direction, 
while in the B layer the spins point opposite to i, as sketched 
in Fig. 1. To avoid confusion with the original bosonic 
operators, the Holstein-Primakoff bosons on layers A and B 
are denoted by a\ , and h\^hi. The resulting Hamiltonian is 
quadratic, and can be diagonalized by a Fourier transforma- 
tion followed by a Bogoliubov transformation, which leads to 
the diagonal Hamiltonian: 

Hlsw = (^k^c^k + ^k/^i/^k) + const., (7) 

k 

where the Bogoliubov operators a\., f3^ are linear combina- 
tions of the Hostein-Primakoff operators aj,ar and 6j,6r- 
The new GS is now the vacuum of the Bogoliubov quasi- 
particles: (o^kO^k) = (/^k/^k) ~ ^' inverting the relation 
between the Bogoliubov operators and the Holstein-Primakoff 
operators gives access to the GS properties of the system. 
The A layer is almost unaffected by the presence of the B 
layer: the density grows linearly above pc (strictly speaking, 
additional small logarithmic corrections are expected in 
2D), as in an isolated layer, and condensate and superfluid 
densities are almost equal. The physical properties of the B 
plane are very different however from those of an isolated 
layer. To describe them more precisely, we consider the 
total density of bosons pf^^ = Xlk(^k^k)' condensate 
density Pq = ;^(^kQ^ko)' superfluid stiffness T^, 

defined by introducing a small twist angle (p along one 
direction of the B plane only. The hopping amplitude t in 
[Eq. (5)] is replaced by te^"^ along this direction, and 
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FIG. 2: Momentum distribution of the LSW modes on the B layer: 3. ^h^^^ ^^^^^ ^j^^ ^-^^^ 



A:^) fort = 1, = 0.5, A = 1, and /i - /ic = A/2 

the superfluid stiffness is given by = [d^ej d^p^^ 
e being the energy density of the system. This leads to the 
expression: = — X]k(^k\) (cos/c^; + cosky). A key 
quantity here is the momentum distribution n^(k) = {bl.b^) 
that we show in Fig. 2 for a representative set of parameters. 
It displays 4 incommensurate maxima away from the BEC 
point ko at Q* = (tt ± g*, tt ± g'*) whose position in the 
Brillouin zone (BZ) depends on A and /i. Perhaps more 
importantly, the momentum distribution vanishes when 
k ^ ko as n^(k) ~ ||k — ko||^ as we discuss below. This 
remarkable behaviour results in the following properties for 
the bosonic fluid in B: (i) the total density pf^^ is non zero 
as soon as /i > /ic [24]; (ii) the condensate mode remains 
empty: quantum fluctuations do not change the classical 
result (^kQ^ko) = ill the GS and the distribution n^(k) 
is never singular, which signals the absence of condensate 
at any vector k; (iii) there is a finite superfluid stiffness 
which is strongly influenced by the location in the BZ of the 
incommensurate vector Q* where the response is maximal. It 
is positive and leads to a superfluid density pf = /2t of 
the same order as the total density if Q* is close to the BEC 
point ko. However, it strongly decreases when Q* shifts away 
from ko, and it eventually changes sign for large enough A 
when the maximum gets far enough from ko, which signals 
an instability towards spontaneous currents, as discussed 
some time ago in the context of dirty superconductors [25]. 

Filtering and frustrated proximity effects — The basic 
physical mechanism behind the absence of a condensate and 
the incommensurate fluctuations in the S-layer is actually 
most simply understood as a filtering mechanism due to 
frustration. To get an intuitive picture of this effect, it is 
useful to look at the single particle dispersion 
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ky = k compared to n'^(k) x 
VA-sik) discussed in the text. The parameters used are t = 1, 
= 0.1, A = 1 at various distances from the critical point /i — /ic- 



ping amplitude t± cos (^) cos (^^^ which vanishes exactly 
at the BEC point ko. Therefore, a single particle injected in 
the A layer with momentum ko is confined to this layer. At 
second order in t± / A, the probability for a boson of momen- 
tum k to tunnel from ^ to >B is given by 



Va-b{^) ^ 



t± ( kx 
— - cos — 
A V 2 
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(9) 



As long as the density in the A layer is not too large, this 
single particle tunneling probability can be used to estimate 
the occupation in S, which, since all particles come from A 
through tunnelling, is expected to be approximately given by: 



^(k) 



(k) x7^^_^(k). 



(10) 



The inter-layer delocalization is governed by an effective hop- 



In the vicinity of the BEC mode ko, the hopping probability 
vanishes very rapidly, as Va-b ^ ll^o — k||^. Combined with 
the fact that the zero-mode fluctuations in the A layer diverge 
like - l/||ko - k||, this leads to n^(k) - ||ko - k|p 
for k ko, in agreement with the exact evaluation of 
n^(k) from the Bogoliubov transformation (see Fig. 3). 
More interestingly, this k-dependent tunneling mechanism 
provides a simple explanation of the maximal response 
observed in B at the incommensurate vector Q*, away 
from the condensate point ko. Indeed, there is a threshold 
vector i<i{ii) {ji — Pc)^^^ such that for ||ko — k|| > the 
occupation in A start to decay faster than l/||ko — k||^, which 
leads to a maximal response at |7r — (5*| ~ The nature 
and the properties of this unconventional quantum liquid are 
summarized in Fig. 4. 

Experimental consequences — Since the filtering mecha- 
nism induced by a frustrated coupling between bosonic layers 
yields a condensate-free superfluid state with incommensurate 
correlations, it is legitimate to ask whether its experimental 
realization is possible or not. As discussed in more details in 
the supplementary material, two directions can be considered, 
namely ultracold bosonic atoms loaded in an optical lattice. 
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qualitatively change the physics. Furthermore, the quantum 
antiferromagnet BaCuSi206 in a field is also predicted to 
display a very rich physics for the triplet excitations in a 
window of ^ 2 Tesla above the critical field He. We hope that 
the present work will encourage experimental investigations 
in both directions, cold atoms and quantum magnets. 
We thank C. Berthier, M. Horvatic and S. Kramer for useful 
discussions on the properties of the Han purple as revealed 
by NMR, and F. Alet and J.-N. Fuchs for a critical reading 
the manuscript. This work has been supported by the Swiss 
National Fund and by MaNEP. 
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FIG. 4: Phase diagram of the bosonic gas of the B layer for the 
frustrated bilayer system, with the normalised chemical potential 
(/i — /ic)/A on the x-axis and the energy barrier A between layers 
A and B on the y-axis. Colours indicate the values of the incommen- 
surate vector Q*. For /j. — /j^c > A = ko = tt, superfluidity and 
EEC occurs together (red region) whereas for /x — /Xc < A, there is 
no true EEC but a finite superfluid density while the stiffness T^, 
positive for not too large A can change sign and become negative for 
larger A (see text). 

and the quantum antiferromagnet BaCuSi206 in a field. 
For the former example, while the harmonic trap will break 
the translational invariance and slightly perturb the perfect 
frustration between layers, we can extend the SW calculation 
to the inhomogeneous case (see supplementary material) to 
show that, in practice, it will not qualitatively change the 
predictions made for the homogeneous case. Concerning 
the frustrated magnet BaCuSi206 in a strong external field, 
as also discussed in the supplementary material, it is a very 
interesting 3D realization of the frustrated layered model 
with two types of layers [26-28]. While the formal analogy 
between superfluids and quantum magnets is limited, we 
predict that high field NMR and neutron experiments may 
detect the absence of field-induced triplets BEC and the 
presence of incommensurate correlations in half the layers of 
the quantum antiferromagnet BaCuSi206. 

Conclusions — We have shown in the context of a simple 
hard-core boson model of two coupled planes that frustration 
can influence dramatically the proximity effect induced by a 
hopping term between the layers. If the frustration is such 
that the condensate of one layer cannot tunnel to an otherwise 
empty layer, then the bosonic gas induced in this layer by 
tunnelling has been shown to be superfluid but uncondensed. 
Beyond this simple model, this effect is expected to be 
present whenever two bosonic systems, an occupied one 
and an empty one, are put into contact by a proximity 
effect, provided the geometry is such that the condensate 
cannot hop from one system to the other. This effect leads 
to non-trivial predictions for trapped bosons in an optical 
lattice where the presence of a quadratic potential does not 
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TRAPPED BOSONS IN OPTICAL LATTICES 



o 
S3 



(N 

a 

o 

(N 
> 
00 

ON 

in 

o\ 
o 
o 



Experiments using ultra-cold gases have recently shown their ability to reproduce strongly correlated many -body physics [1]. 
In particular, when loaded in a D-dimensional optical lattice, trapped bosons may display quantum phase transitions such as the 
celebrated BEC-Superfluid to Mott-Insulator transition [2]. Another remarkable achievement concerns the interesting Tonks- 
Girardeau regime analog to hard-core bosons [3]. Recently, it has also been shown that artificial gauge fields [4, 5] can be created 
in such systems, allowing to achieve positive hopping terms between nearest neighbors — texp(i7r), a necessary condition to 
frustrate the condensate mode in our case. As compared to archetypal solid-state uniform bulk systems, the price to pay with 
optical lattices is double: finite system size (typically 10^ lattice sites) and a trapping potential which breaks the translational 
invariance. The presence of a harmonic trap leads to a site-dependent chemical potential 
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Usually, a method known as local density approximation is used to deal with inhomogeneous potentials [6]. However, this is not 
so simple here since frustration turns out to be quite sensitive to the inhomogeneities in the local densities. Indeed, contrary to 
the translationally invariant case where in the magnetic language a uniform antiferromagnetic XY order on the A layer yields 
a zero local mean field at each site of the B layer, the classical spin vectors S = 5'(sin 6^, 0, cos 0) have angles which are now 
smoothly site-dependent ^^(r) within the trap. As a consequence, ^^(r) (pl.=plaquette) will not necessarily vanish, as 
illustrated in Fig. 1 (d) in a ID representation. 

The classical solution is obtained by a global numerical minimization of the full set of local angles 6j[/]s{t) on both layers. 
Unlike the uniform case (Vb = 0)» here an extremely small but finite population can hop from AtoB already at the classical level, 
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FIG. 1: Real space density profiles for a trapped bilayer system of 24 x 24 x 2 lattice sites, with a chemical potential /i = iiic + t/2,3. trapping 
potential Vb = 10~^t, a barrier A = 4t, and an interlayer hopping t± = 0.5t. Local densities at the classical approximation shown for layers 
A (a) and B (b). Quantum fluctuations induce a much larger population in B (c). (d) Top: schematic picture for both layers A and H in a 
harmonic trap and separated by a barrier A. Down: ID representation of the frustrated setup which gradually loose perfect frustration away 
from the trap center. 
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as visible in Fig. 1 (b). Nevertheless, this classical population has a peculiar spatial distribution due to the trade-off between two 
competing effects. On one hand, in order to satisfy the local energetic constraints imposed by the chemical potential, particles 
would tend to be maximally present at the center of the trap with a density which decreases away for the trap center, as clearly 
observed for the A layer in Fig. 1 (a). On the other hand, the local constraint imposed by the interlay er frustration decreases 
with r. At the very center, and along the axes of the square lattice, the symmetry imposes a perfect frustration ^a{'^) = 
which gradually disappears away from these lines. These two competing effects give rise to a distribution in the B layer which is 
localized away from the center, the square lattice axes, and the edge of the trap, as shown in Fig. 1 (b). Remarkably, the maximal 
density in S is 5 orders of magnitude below the maximal density in the A layer. This shows that the frustration effect, while 
formally unperfect, remains in practice very efficient, especially close to the center of the trap. Nevertheless the nature of this 
highly diluted bosonic population which has tunneled from ^ to S is conventional, i.e. condensed and superfluid but with a tiny 
density, almost experimentally undetectable. 

As previously done for the uniform case, we investigate quantum fluctuations on top of the classical GS using the linear 
spin- wave approximation. Due to the lack of translational invariance, we must study quantum fluctuations in real space, which 
leads to an eigenvalue problem of much larger dimension (twice the number of sites instead of 4). The GS energy Eo{{jj.{r)}) 
for a given set of chemical potentials is obtained as the sum of the classical energy and of the zero point fluctuations, and the 
local density at a site r is obtained computing —dE^/diiiY) evaluated numerically at the chemical potentials of Eq. (1). As 
visible in Fig. 1 (c), the bosonic fluid induced by quantum fluctuations completely dominates the physics of the B layer, in which 
the classically empty center of the trap is now maximally populated. Most importantly, this new bosonic population, 3 order 
of magnitude above the classical one, will display similar properties as in the uniform case, i.e. absence of Bose condensate, a 
finite superfluid density, and incommensurate correlations. Newly developed detection methods in real space [7, 8] will allow 
to probe the bosonic population close to the trap center. The incommensurate response at Q* can be detected by measuring the 
momentum distribution n^(k) using standard time of flight experiments [3]. 



BACUSI2O6 AT HIGH MAGNETIC FIELD 

As already discussed in Ref. [9-11], BaCuSi206 under a strong magnetic field is a good realization of a 3D arrangement of 
inequivalent bosonic planes A-B which display a fully frustrated transverse coupling, and a bare gap A 3 T. Building on our 
understanding of the two-layer boson system studied above, it is straightforward to explore the properties of the exotic triplet 
fluid on B layers within a LSW approach of the full 3D model with realistic parameters for BaCuSi206. In a magnetic field 
range of ~ 2 T above the critical field He = 23 A T, the B layers will not contribute a Bragg peak at ko = (tt, tt), but rather 
develop a maximum at a wave- vector Q* which is expected to display an unusual power-law dependence on the distance to the 
critical field He with an exponent 1/4: 

\\Q*-U\^{H-H,)'/\ (2) 

These features are direct consequences of the absence of condensate, one of the main predictions of the theory developed in this 
paper. They are in principle accessible to NMR and neutron scattering. In NMR, the absence of condensate implies that the 
massive broadening of the line that accompanies the development of magnetic order will not take place. The incommensurate 
response can in principle be observed in the line shape. This will be complicated by the presence of an incommensurate distortion 
in the Han purple [12], but the two effects might be disentangled by studying the evolution of the line shape with the field since the 
incommensurate response predicted by the present theory is field dependent, while the incommensurate distortion is not expected 
to be. The absence of a Bragg peak and the incommensurate response may also be detected by neutron scattering. Indeed, using 
pulsed field, the possibility to perform elastic neutron scattering up to 30 T has been very recently demonstrated [13], thus 
opening the route to probe the peculiar field-induced momentum distribution in BaCuSi206. Let us reemphasize however that 
the other prediction of our theory, the persistence of superfluidity despite the loss of condensate does not apply to the Han purple 
because of residual anisotropics [14] that break the rotational symmetry. 



[1] L Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008). 

[2] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and L Bloch, Nature 415, 39 (2002). 

[3] B. Paredes, A. Widera, V. Murg, O. Mandel, S. Foiling, L Cirac, G. V. Shlyapnikov, T. W. Hansch, and L Bloch, Nature 429, 277 (2004). 

[4] Y. Lin, R. L. Compton, K. Jimenez-Garcia, J. V. Porto, and I. B. Spielman, Nature 462, 628 (2009). 

[5] J. Dalibard, F. Gerbier, G. Juzeliunas, and P Ohberg, ArXiv e-prints (2010), 1008.5378. 

[6] G. G. Batrouni, V. Rousseau, R. T. Scalettar, M. Rigol, A. Muramatsu, P J. Denteneer, and M. Troyer, Phys. Rev. Lett. 89, 1 17203 (2002). 



3 



[7] J. R Sherson, C. Weitenberg, M. Endres, M. Cheneau, L Bloch, and S. Kuhr, Nature 467, 68 (2010). 

[8] W. S. Bakr, A. Peng, M. E. Tai, R. Ma, J. Simon, J. L Gillen, S. Foiling, L. Pollet, and M. Greiner, Science 329, 547 (2010). 

[9] C. Riiegg, D. F. McMorrow, B. Normand, H. M. R0nnow, S. E. Sebastian, 1. R. Fisher, C. D. Batista, S. N. Gvasaliya, C. Niedermayer, 
and J. Stahn, Phys. Rev. Lett. 98, 017202 (2007). 
[10] S. Kramer, R. Stern, M. Horvatic, C. Berthier, T. Kimura, and I. R. Fisher, Phys. Rev. B 76, 100406 (2007). 
[1 1] N. Laflorencie and F. Mila, Phys. Rev. Lett. 102 (2009). 

[12] E. C. Samulon, Z. Islam, S. E. Sebastian, P. B. Brooks, M. K. McCourt, Jr., J. Ilavsky, and L R. Fisher, Phys. Rev. B 73, 100407 (2006). 
[13] S. Yoshii, K. Ohoyama, K. Kurosawa, H. Nojiri, M. Matsuda, P. Frings, F. Due, B. VignoUe, G. L. J. A. Rikken, L. Regnault, et al., Phys. 
Rev. Lett. 103, 077203 (2009). 

[14] S. E. Sebastian, R Tanedo, R A. Goddard, S. Lee, A. Wilson, S. Kim, S. Cox, R. D. McDonald, S. Hill, N. Harrison, et al., Phys. Rev. B 
74, 180401 (2006). 



